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Abstract 

We discuss the implementations of the Bose-Einstein effect from asymmet- 
ric sources in Monte Carlo generators. A comparison of LEP data with results 
from the PYTHIA/JETSET code with the standard procedure imitating the 
effect and with the results from the weight method (with weights depending in 
various ways on components of momenta differences) is presented. We show 
that in this last method one can reproduce the experimental hierarchy of the 
source radii. 



1 Introductory remarks 



Recently one observes a renewal of interest in analysing the space-time structure 
of sources in multiparticle production by means of Bose-Einstein (BE) interference 
. Such analysis followed the example of astrophysical investigations of Hanbury- 
Brown and Twiss . The main motivation of this renewal was the analysis of the 
e + e~ — > W + W~ process which became available at the LEP2. It was suggested 
0,0 that the BE interference (and/or colour reconnection effects) between the 
strings from two W decays may shift the W mass value fitted from the two jet mass 
distributions by as much as a few hundred MeV, thus making this channel useless 
for precise tests of the standard model. However, other investigations suggested that 
such a big shift is unlikely 0,0,0. Experimentally, the existence of interference 
effects between strings is still debatable [§]. 

Investigating such subtle effects became possible when instead of the standard 
approach one started to model this effect in Monte Carlo generators. There 
are several methods of modelling: as the "afterburner" for which the original MC 
provides a source [10|,[11]], by shifting the momenta Jl2 or by adding weights to 
generated events [13],[14|. Another approach was set forward by Andersson and 



collaborators who used the symmetrization inside fragmenting string [[Uj to model 
the effect for a single string (T^|. Here we consider the most widely used methods of 
shifting momenta and weighting events. 

Another reason to analyze the BE effect were the efforts to estimate size and 
shape of source of particle production in various processes (in particular for coming 
RHIC data). The analysis of BE effect in 3 dimensions is supposed to reflect the 
spatial source asymmetry. Such analysis was done for the LEP data at the Z° 
peak |17j which have very high statistics and good accuracy. 

In this paper we compare the 3-dimensional data for BE effect from LEP with 
the results of the standard momentum shifting procedure and of the weight method. 
In the next section we present the data discussing in detail the definitions and 
the procedures used by the experimental groups. In the third section we compare 
them with the results obtained from the PYTHIA/JETSET MC generator using the 
original procedure modelling the effect by momentum shifting and with the results 
from the weight method with weights independent on spatial orientation of momenta. 
Fourth section contains the results for asymmetric weights. Our conclusions are 
presented in the last section. 



2 Experimental data 

Although the discussion of the shape of asymmetric sources in the framework of BE 
interference concerned most often the heavy ion collisions, the best experimental 
data with highest statistics exist for the e + e~ annihilation at the Z° peak. In the 
following we concentrate our attention on the L3 data [Tj| which discuss the ratios 



using "uncorrelated background" and three different radii to parametrize the data. 
The DELPHI data |B1 are parametrized with only two radii, and the OPAL data 
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use the like/unlike ratio which requires a cut off of the resonance affected regions 
even in double ratios. 

As in the L3 paper jLl| we use for each pair of identical pions three components of 
the invariant Q 2 = — (pi — P2) 2 '- Q 2 L,Qout,Q 2 ide defined in the LCMS (Longitudinal 
Center-of-Mass System), where the sum of three - vector momenta is perpendicular 
to the thrust axis. The Q out component is measured along this sum, the Ql along 
the thrust axis, and Q S ide is the projection of Q on the axis perpendicular to these 
two directions 0,0. 

We define a "double ratio" in the same way as in the L3 paper using a reference 
sample from mixed events: 

P2 P2 IC 

R 2(Pl,P2) = -—/ mixM C- W 
P2 P2 

This "double ratio" is parametrized by 

R2^Ql, Qout, Qside) = 7[1 + &Ql + ^Qout + CQ side] ' 

■ [1 + A exp(-R 2 L Q 2 L - R 2 out Q 2 out - R 2 side Q 2 side - 2pL,outRLRoutQLQout)}- (2) 

The first bracket reflects possible traces of long-distance correlations; the last 
term in the second bracket seems to be negligible when fitting data and will be 
omitted in the following. 

By fitting the parameters Rl and R S id e we get some information on the geometric 
radii in the longitudinal and transverse directions (respective to the thrust axis). 
R out reflects both the spatial extension and time duration of the emission process. 

In the L3 data the fit region in all three variables extends to 1.04 GeV and 
is divided into 13 bins, which gives 2197 points fitted with 8 parameters. The fit 
parameters 5, e and ( are rather small; this means that the observed BE enhancement 
is rather well approximated with a Gaussian. The value of the parameter A is fitted 
as 0.41 ±0.01. 

The fitted values of radii (in fm) are as follows: 

R L = 0.74 ± 0.02jft{j|, R out = 0.53 ± 0.02jft{jjj, R side = 0.59 ± 0Mt° H 

We see clear evidence for source elongation: R S idelRh is smaller than one by more 
than four standard deviations. 

It is instructive to inspect the projections of the double ratio on the three axes 
Ql, Qout and Q S ide- This is done by restricting the values of two other variables to 
less than 0.24 GeV, plotting the histograms in the third variable in bins of width 
0.08 GeV and constructing the double ratio in this variable. The results are shown 
in Fig. |1] as presented by the L3 collaboration 



The values of double ratios fall down smoothly from the maxima of about 1.25 
at Q-s close to zero to the plateau at 1. It is rather difficult to see the differences 
between three plots, but superposing them one may note that the fall is fastest for 
Ql as expected from the fact that the fitted value of parameter R L is bigger than 
the values of R out and R S ide quoted above. 
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Figure 1: Projections of the double ratio (|IJ) from the data of the L3 collaboration 
on the three axes Ql, Q ut and Qside- 

3 Asymmetric effects from symmetric models 

The geometric interpretation of data requires a comparison with the results from the 
standard MC procedures modelling the BE effect. In the L3 paper such an analysis 
is given for the standard LUBOEI procedure built into the JETSET Monte Carlo 
generator. This procedure modifies the final state by a shift of momenta for each 
pair of identical pions. The shift is calculated to enhance low values of Q 2 and to 
reproduce the experimental ratio in this variable. The function defining this shift is 

f(Q 2 ) = l + X m exp(-R 2 n Q 2 ) (3) 

The superposition of the procedure for all the pairs and subsequent rescaling 
(restoring the energy conservation) makes the connection between the parameters 
of the shift Aj n , Ri n and the parameters describing the resulting double ratio in Q 2 

Pi p¥° 

Rl{Q ) = m ix,MC (4) 

P2 p 2 

(which may be parametrized analogously to (|3|)) rather indirect. 

Using the JETSET parameters adjusted to all the L3 data and the LUBOEI 
parameters fitted to describe the BE ratio in Q 2 the authors of the L3 paper cal- 
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culated the same quantities as measured in the experiment. The projections of R2 
are qualitatively very similar to the experimental ones. However, the fit to the 3- 
dimensional distribution gives results different from data. The ratio R S id e /RL is not 
smaller but greater than one; the fitted values (in fm) are: 

R L = 0.71 ± 0.01, R out = 0.58 ± 0.01, R side = 0.75 ± 0.01. 

We confirmed these numbers in our calculations. We found also that the results are 
sensitive to the JETSET parameters. Using the default values instead of the L3 
values we obtained a significantly smaller value of R out (below 0.5) and significantly 
smaller A. Other values are less affected and R S id e /Ri is still bigger than 1. 

We have also checked how the results depend on the source radius Ri n and on the 
incoherence parameter X in assumed in the LUBOEI input function @. In all cases 
we get R S ide > Rl > Rout, although the input function was obviously symmetric. 
The values of R S i de and Rl are proportional to R in , whereas R out changes much less; 
the dependence on X in is very weak. The output value of A decreases quite strongly 
with increasing Ri n and increases with Aj„. No choice of input parameters gives the 
values of Ri compatible with data. This is shown in Fig. |^. 
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Figure 2: Fit parameters A and Ri as functions of the input parameters a) for the 
LUBOEI procedure, b) for the weight method. Experimental values are shown on 
a separate vertical axis. 

Another interesting observation is that to fit the L3 data one needs A = 1.5, 
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which is beyond the physically acceptable value of 1. This supports our doubts 
about usefulness of the LUBOEI procedure in understanding the experimental re- 
sults (although certainly it is the most practical description of data). 

In fact, there is one more degree of freedom in the prescription for modelling 
the BE effect: the definition of direct pions. Since the decay products of long-living 
resonances and of particles decaying by electroweak interactions are born far from 
the original collision point, their effective source size is much bigger than that for 
direct pions. Thus they contribute to the BE effect for momentum differences much 
below the experimental resolution and should not be taken into account. 

In the LUBOEI procedure this distinction is made by the decay width of unstable 
particles: only pions from the decay of particles with the width above 20MeV and 
the direct ones are included in the momentum shifting procedure. Obviously, this 
is just a rough prescription which may be changed, and the values of fit parameters 
may change then quite strongly. The user of the procedure should be aware that 
(according to author's warning) it works properly only when called from LUEXEC; 
if LUBOEI is called directly from the master program, all pions are regarded as the 
direct ones. 

The problems of LUBOEI procedure in describing the asymmetry of experimental 
distributions are not the first ones noted in applications to describe various data. 
It has been already indicated that the procedure with parameters fitted to the two- 



particle data fails to reproduce the three-particle spectra |22| and the semi-inclusive 



data [23]. Moreover, as already noted, the fitted values of parameters needed in 
the input function (H) are quite different from the values one would get fitting the 
resulting double ratio (f|) to the same form p3J , PI| . Thus, it seems to be difficult 



to learn something reliable on the space-time structure of the source from the values 
of the fit parameters in this procedure. 

All this led to a revival of weight methods, known for quite a long time |26 



but plagued also with many practical problems. The method is clearly justified 
with in the formalism of the Wigner functions, which allows one to represent (after 
some simplifying assumptions) any distribution with the BE effect built in as a 
product of the original distribution and the weight factor, depending on the final 



state momenta [13]. With an extra assumption on factorization in momentum space 



we may write the weight factor for a final state with n identical bosons as 

n 

W(pi,...p n ) =J2Y[ W 2(Pi,PP(i)), (5) 

i=l 

where the sum extends over all permutations P n {i) of n elements, and W2{pi,Pk) is 
a two-particle weight factor reflecting the effective source size. Problems with an 
enormous number of possible terms in this sum may be cured by a proper clustering 
procedure WM. A reasonable description of the effect in Q 2 is obtained with a simple 



gaussian form of the weight factor 

w 2 (pi,p 2 ) = exp[-(p! -p 2 ) 2 R 2 n /2], (6) 
or, even simpler, a step function form with w 2 = 1 for some range of —(pi — p 2 ) 2 < 



l/Rf n and w 2 = outside |27 
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In this method we may repeat the same calculation as done for the LUBOEI 
procedure. Obviously the weights may be calculated for the events generated by 
any MC generator, but here we restrict ourselves to the results from the same 
PYTHIA/JETSET code which was used above. The resulting double ratios are 
not that smooth and monotically decreasing as in the data or from the LUBOEI 
procedure (which is the usual drawback of the weight methods). However, the major 
features are surprisingly similar: with weight factors depending only on Q 2 we get 
different values of fitted Ri parameters. Moreover, the hierarchy of parameters is 
the same: R S id e > Rl > Rout- This suggests that the assymetry is generated by the 
jet-like structure of final states and not by any specific features of the procedure 
modelling the BE effect. In Fig. 2b we show the values of the fit parameters as 
functions of Ri n for a Gaussian as well as the #-like weight factors. Again, no choice 
of the input parameters allows to describe the data. 

The comparison of two methods is not straightforward. In particular, one should 
take care if the same definition of "direct" pions is used. The weights are calculated 
after the event was fully generated (and all the decays of unstable particles occured). 
Therefore one should define the pions which are counted as direct ones. We did it by 
enumerating particles, which contribute significantly to the pion production and live 
too long for their decay products to produce a visible BE effect (using the same limit 
for decay width as in LUBOEI). If one enumerates the short-living resonances and 
adds their decay products to the direct pions, one should remember that this list is 
different in various options of JETSET (e.g. the option used by the L3 collaboration 
takes into account mesons built from quarks with non-zero orbital momentum, which 
are neglected in the default version). 

The results presented in this section suggest that one should be careful with the 
geometric interpretation of the data. If one gets asymmetric distributions from the 
generator without assuming explicitly space asymmetry of the source, it is not clear 
how the assumed asymmetry will be reflected in the results. 

4 Asymmetric weights 

One may get more information on the problem of asymmetric BE effect in MC 
generators using the asymmetric weight method, i.e. introducing weight factors 
which depend in a different way on Q L = \p 1L - p 2 i\, Qside = \puide - P2side\ and 
Qout = \p\out— Viout I, where the indices denote the components defined in the previous 
section. We have used two such generalizations of a gaussian weight factor (|]) 

w 2 (Q L , Q^, Q side ) = exp([-gi( J R™) 2 - Q^Kutf - QLeiKlde) 2 } 1 2) (7) 
and 

MQl, Qout, Qside) = exp([-Ql(i?l") 2 - (1 - (3 2 )Q 2 ou Mlt) 2 - QLefe) 2 ]/2) (8) 
where (3 is defined as 
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n _ VoutX + Pout2 /„n 

P ~ E 1 + E 2 ■ (9) 

The weight factor (f|) reduces to the symmetric weight factor (|6|) when R l £ = 
R % out — R % Iide = Rin- The formula @ gives nearly the same results as the formula (0) 
when R out is multiplied by 2. We have used both forms finding no definite preference 
for any of them. 

Fluctuations in the weight values are large and the resulting fluctuations in the 
values of double ratios describing the BE effect are bigger than for the momentum 
shifting method. Therefore it is necessary to use large samples of generated events. 
We found that for the samples of 5 million events the fluctuations visible in the 
plots of projections of double ratios on components of Q are comparable with those 
seen in the experimental data shown in Fig.l. In fact, the plots obtained for the 
weight method with the input radii around 0.5 fin are visually similar to those of 
experimental data. However, the fitted values of the parameters from formula @ 
are different. 

Since for the symmetric weights the resulting fitted values of R S ide are bigger than 
the values of Rl (contrary to the inequality seen in the data), it seemed natural to 
take the input value of R™ de smaller than R % £. Indeed decreasing R™ de one reduces 
the resulting fitted value of R S id e , but this dependence is not linear and saturates 
for R™ de around 0.3 fm. Moreover, the fitted values of other parameters change as 
well although their input values were not changed. Therefore looking for the best 
set of input parameters in the formula for weights is a rather involved procedure. 

Let us add two more remarks. A replacement of the products of Gaussians by 
the proper products of step functions in the formulae for weights (|7p, ([5p leads to 
even bigger fluctuations in the resulting distributions and we do not advocate such 
parametrizations. Finally, there is some ambiguity concerning the use of weights for 
the calculations of double ratio ([[]). If we use the weights only for the two-particle 
distributions, the two denominators cancel and we calculate effectively just the ratio 
of two-particle distributions with- and without weights. It seems, however, that the 
justification for the weight method |TB[ requires using weights both for the single- 



and two-particle distributions. We have looked for the best set of parameters with 
this prescription, using a Gaussian form without the "/5-factor" (|7|). The best set 
we found is 

R? = 0.9 fm, RZ t = 0.3 fm, Rf ide = 0.4/m. (10) 

The resulting projections of the double ratios are shown in Fig. 3. The fitted 
values of parameters in formula (0) we get are 

R L = 0.73 fm, R out = 0.54/m, R slde = 0.65/m. (11) 

Obviously, it is now possible to reproduce the experimental hierarchy of radii. 
The fitted value of A is smaller than in data (0.35 instead of 0.41), but the difference 
is well within the systematic errors of the fit to the experimental data. Note that we 
are not showing the errors in Fig. 3 (nor quoting them in the values of parameters 
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listed above), since these errors result mainly from the fluctuations in weights. Some 
estimate is obtained by comparing the results for 1 and 5 million events samples; in 
Fig. 3 the differences are of the order of size of the points. 
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Figure 3: Projections of the double ratio (g) from the PYTHIA/JETSET MC gen- 
erator with the asymmetric weight method for parameters (|TT|) on the three axes 
Qli Qout and Q side- 
There is a striking difference between the input values of the radii fllOD assumed 
in the weight factors and the resulting best fit values (|TTJ) from the double ratio 
calculated with these weights. Although the hierarchy Rl > R S id e > Rout is the 
same in both cases, the fitted values differ by less than 25%, whereas there is a 
difference by more than a factor of two between the input values. 

Moreover, further decrease of the values of R^ t and R l ^ de hardly affects the 
resulting double ratio and fitted values of R4. This seems to be the inherent property 



of the JETSET generator, which yields a rather strong suppression of large values 
of Qi and Q 2 even without any procedure imitating the BE effect. Apparently 
this suppression dominates over the weak enhancement of low values of Qi induced 
by the weight factors with small values of R{. For small _R- n there is no simple 
correspondence between the input and output values of radii. This looks analoguos 
to the effect noted already for a symmetric BE effect described by the LUBOEI 



procedure p5| . Therefore any direct interpretation of the fit values for BE double 
ratios in terms of the different radii of the asymmetric source is a rather delicate 
matter. 



5 Conclusions 

In this note we present the results of our investigation concerning the asymmetry 
of the BE effect in two procedures imitating this effect in the Monte Carlo genera- 
tors. A comparison with the data at Z° peak is presented. We found that both the 
momentum shifting method and the weight method with weights depending on Q 2 
only give different distributions in different components of Q 2 . However, the hierar- 
chy of radii parametrizing these distributions is different from the experimental one. 
Introducing weights which depend in different ways on different components of Q 2 
we are able to reproduce the experimental data. 
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